Anderson localization of matter waves in tailored disordered potentials 
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We show that, in contrast to immediate intuition, Anderson localization of noninteracting particles 
induced by a disordered potential in free space can increase (i.e. the localization length can decrease) 
when the particle energy increases, for appropriately tailored disorder correlations. We predict the 
effect in one, two and three dimensions, and propose a simple method to observe it using ultracold 
atoms placed in optical disorder. The increase of localization with the particle energy can serve to 
discriminate quantum versus classical localization. 
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The transport properties of a coherent wave in a disor- 
dered medium are inherently determined by interference 
of multiple scattering paths, which can lead to spatial lo- 
calization and absence of diffusion [1]. This effect, known 
as Anderson localization (AL), was first predicted for 
electrons in disordered crystals [2] and then extended to 
classical waves [3] , which permitted observation of AL in 
a variety of systems (sec Rcf. [4] and references therein). 
The most fundamental features of AL are ubiquity and 
universality [5]. For instance, in conventional cases, all 
states are known to be localized in one (ID) and two (2D) 
dimensions, while in three dimensions (3D) the spectrum 
splits into regions of localized states and regions of ex- 
tended states, separated by so-called mobility edges [6]. 
Nevertheless, the observable features of AL strongly de- 
pend on the system details. 

Consider a wave propagating among randomly dis- 
tributed point scatterers (point-impurity disorder). In 
the absence of interference, the propagation is dominated 
by normal diffusion. It devises a diffusive medium char- 
acterized by the length scale (transport Boltzmann mean 
free path) l B = vt, with v = \du(k)/dk\ the wave veloc- 
ity [w(fc) is the dispersion relation] and t the scattering 
time. Then, localization arises from the interference of 
the diffusive paths. The more the wavelength exceeds 
the mean free path, the stronger interference affects the 
transport. It can thus be inferred that the Lyapunov ex- 
ponent (inverse localization length), which characterizes 
the localization strength, reads 7 = i" 1 Fd{kl B ) where the 
function strongly depends on the spatial dimension d 
and is a decreasing function of the interference parameter 
kl B . For a particle (scalar matter wave) in free space and 
a weak point-impurity disorder, v cx k, t is proportional 
to the inverse of the density of states (p cx k d ~ 2 ), as given 
by the Fermi golden rule, and finally l B oc l/k d ~ 3 . Then, 
for any d < 3, 7 is a decreasing function of k. In other 
words, the localization gets weaker when the particle en- 
ergy E = huj(k) = h 2 k 2 /2m increases, which conforms to 
natural intuition. 

This decrease of 7(E) with E however relies on the mi- 
croscopic details of the system, namely on the dispersion 
relation and on the properties of the scattering time as- 
sumed above, and can be altered in different ways. For 



instance, it does not hold for lattice systems, such as elec- 
trons in disordered crystals, because the band structure 
leads to a nonmonotonic behavior of v versus E, which 
can lead to a function 'y(E) approximately symmetric 
with respect to the band center [7, 8]. In other systems, 
such as light waves in dielectric media, r shows Mie res- 
onances [9], leading to a strongly nonmonotonic behav- 
ior of j(E). In this work, we discuss a different effect. 
We show that the standard behavior of j(E) for parti- 
cles in free space can be inverted (i.e. localization can get 
stronger with increasing energy) by tailoring the disorder 
correlations. The basic idea behind our work is that for 
non-point scatterers, the structure factor C{k) appears 
in the denominator of r. Then, if the disorder has strong 
spatial frequency components around a particular value 
fco, the scattering strength may not vary monotonously 
with E around E(ko), and j(E) can then increase with 
E. In contrast to the cases discussed above (lattice elec- 
trons and light waves in dielectric media), this effect is 
purely due to the disorder correlations. We first study 
the ID case, which allows for exact calculations of 7(2?) 
and for explicit test of an efficient scheme to observe the 
effect with ultracold atoms. We then extend our analysis 
to 2D and 3D systems using the self-consistent theory of 
AL. We finally discuss how the increase of 'y(E) with E 
can serve to discriminate quantum versus classical local- 
ization. 

A key ingredient of our work is the possibility of tai- 
loring the disorder correlations. Let us consider speckle 
potentials [10-14]. A speckle field [15] is obtained from 
a coherent laser beam diffracted by a ground-glass plate 
diffuscr, which imprints a spatially random phase on the 
electric field at each point R of its surface [see Fig. 1(a)]. 
The diffracted complex electric field £(r) at a given ob- 
servation point r is the sum of independent random vari- 
ables, corresponding to the components originating from 
every point R of the plate and interfering in r. The atoms 
are subjected to a potential, which, up to an arbitrary 
shift, is proportional to the light intensity X cx \£\ 2 - We 
define V(r) = V R x {l(r)/T - 1}, so that V = and 
V 2 = V 2 (note that I 2 = 21 2 [15]). The sign of V R 
can be positive or negative depending on the detuning of 
the laser with respect to the atomic resonance. In the 
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Figure 1. (Color online) Tailoring correlations in speckle potentials, (a) Optical apparatus: A laser beam is diffracted by a 
ground-glass plate diffuser (D) of pupil function 7 D (R), where R = (Y, Z) spans the diffuser, which imprints a random phase 
on the various light paths. The intensity field, T(r), observed in the focal plane of a converging lens, is a speckle pattern, 
which creates a disordered potential V(r) for the atoms, (b) One-dimensional speckle potentials realized with a pupil function 
obtained with two incident Gaussian beams of waist w and centered at Z — ±A/2. The figure shows a sketch of I D (1 st 
column), a realization of V(z) (2 nd column), and the reduced disorder spectrum in fe-space, c(ka R ) (3 rd column) for various 
values of A/w. 



paraxial approximation for the scheme of Fig. 1(a), the 
disorder power spectrum (Fourier transform of the disor- 
der correlation function [16]) in the focal plane (y,z) of 
the lens is [15] 

C*(k)cx J dR / D [R-(A Z/47r)k]/ D [R+(A Z/47r)k] (1) 

where I D (R) is the pupil function (i.e. the intensity pro- 
file right after the diffusive plate), A is the laser wave- 
length, and I is the focal length. The major constraints 
on C(k) follow from Eq. (1) and from the fact that 7 D (R) 
is nonnegative and of finite integral. Firstly, the Cauchy- 
Schwarz inequality applied to Eq. (1) shows that C(k) 
is a decreasing function of |k| for small values of |k|. 
Secondly, in practice, I D (R) decays at long distance so 
that C(k) also decays in the large |k| limit. Apart from 
these constraints, control of / D (R) offers freedom for tai- 
loring the disorder power spectrum, which we will write 
C(k) = V 2 a R c(ka R ) with a R the correlation length. We 
now show that it allows us to strongly affect the qualita- 
tive behavior of AL for noninteracting particles. 

To start with, consider the ID case. Exact calculations 
can be performed [17] within the Born approximation, 
valid for weak disorder [i.e. for 7(E) <C fc^er^ 1 where 
k E = \/2mE /h] . They yield the Lyapunov exponent 

<y(E) = L-j(E) ~ {m 2 V 2 a R /2h A kl)~c{2k E a R ). (2) 

As can be explicitly seen in Eq. (2), if the disorder fea- 
tures no particular correlations, i.e. if c(k) is a constant 
or decreasing function of k, then 7(E) decreases mono- 
tonically with E, and the localization is weaker for higher 
energy, as for point-impurity disorder. In order to invert 
this behavior in a given energy window, it is necessary to 
tailor the disorder correlations so that c(2k E a R ) increases 
with k E strongly enough to overcome the 1/k 2 decrease 
of the prefactor in Eq. (2). To do so, we propose to use 



speckle potentials realized by illuminating the diffusive 
plate by two mutually coherent Gaussian laser beams of 
waist w along Z and centered at Z = ±A/2 [18]. Using 
Eq. (1), we find 



c(k) 
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(3) 



with cr R = \ 1/ttw and k = 2 A/w, the values of which 
can be independently controlled. The properties of the 
disordered potentials obtained in this configuration are 
shown in Fig. 1(b) for various values of k„. For k = 
(lower row), the disordered potential features structures 
of typical width cr R in real space (central column). The 
corresponding power spectrum C(k) has a single Gaus- 
sian peak of rms width y/2/a n centered in k — (right 
column). For k ^ 0, the disordered potential develops 
additional structures of typical width a B /n„, correspond- 
ing in C(k) to an additional peak centered in k ~ K /a R . 
For k large enough, C(k) shows a clear increase with 
A; in a significant range (upper rows). For k > 5.3, we 
find that it is strong enough that -f(E) is nonmonotonic, 
hence realizing the desired situation where localization 
becomes stronger when the particle energy increases. For 
instance, for n — 8.88, 7(E) shows a significant increase 
between k E ~ 2.3CTR 1 and k E ~ 4.2CT" 1 [see Fig. 2(a)]. 

We now discuss how to observe this effect using ul- 
tracold gases expanding in the disordered potential. In 
the long-time limit, for vanishing interactions and an ini- 
tial gas of negligible width, the average spatial density 
reads [19, 20] n(z,t-> 00) = J dE V E (E)P 00 (z\E), where 
T> E (E) is the energy distribution of the atoms and 

r 1 -l ,,2 -> 2 

du u sinh(7ru) 
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xexp{-(l + U 2 ) 7 |z|/2}, (4) 

with 7 = 7(E) given by Eq. (2), is the probability of 
quantum diffusion [21]. Using the scheme of Ref. [10], 
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Figure 2. (Color online) Anderson localization in ID speckle 
potentials with the autocorrelation function (3), k = 8.88 
and Vr = ±0.72(ft 2 /m<r R ). (a) Lyapunov exponent versus 
particle energy as obtained from Eq. (2) (solid black line) 
and from fits of Eq. (4) to numerical data (points), (b-c) Sta- 
tionary density profiles obtained numerically using the initial 
state (5) with p w <j R //i = 0.24 and two different values of p at - 
The figures show the results for six realizations of the dis- 
order [three with Vr > (blue upper data) and three with 
Vr < (red lower data)], the averaged density profile (black 
central data) and the fits of Poo (z) to the latter (green solid 
line). The extracted values of 7f it for each realization and for 
the averaged profile are reported in (a). The figure shows a 
significant increase of the Lyapunov exponent in the energy 
window such that 2.3 < fc B a" R < 4.2. 



for which the energy distribution extends from E = 
to E = i? max , does not allow us to probe the region 
where ^(E) increases because the long distance behavior 
of n co (z) would always be dominated by the energy com- 
ponents with the largest localization lengths, i.e. those 
with the smallest ^(E) [19, 20]. Instead, we propose to 
use an atomic energy distribution strongly peaked at a 
given energy E at , so that n OQ (z) ~ Poo{z\ E &t)- It can 
be realized by either giving a momentum kick to a non- 
interacting initially trapped gas or using an atom laser, 
both with a narrow energy width. The momentum dis- 
tribution can be represented by a ID Gaussian function 
of width p w centered around a controllable value p at [22] : 



V p (p) = (1/V27rp w ) exp \-{p - Pat y/2pi 



(5) 



For weak disorder, the corresponding energy distribution 
is weakly affected by the disorder-induced spectral broad- 
ening, so that it is strongly peaked at E at ~ plJ2m [20]. 

We have performed numerical integration of the time- 
dependent Schrodinger equation for a particle in the dis- 
ordered potential with the initial momentum distribu- 
tion (5) and disorder parameters as in Fig. 2(a). Dur- 
ing the expansion, back and forth scattering processes 
quickly redistribute left- and right-moving atoms. The 
center of the cloud hardly moves and the wings gradu- 
ally form a nearly symmetrical stationary density pro- 
file rioo(z), shown in Figs. 2(b) and (c) for two values of 



p at and for six realizations of the disordered potential: 
three with blue detuning (V R > 0) and three with red 
detuning (V R < 0). The density profile averaged over 
the six realizations, n^(z), is also displayed (black line). 
After averaging, we fit lnfPo^z)] as given by Eq. (4) 
to ln[noo(2:)] with 7 as the only fitting parameter. Al- 
though the fits are performed in a limited space window 
(— 300fi R < z < +300er R , corresponding to an experimen- 
tally accessible width of 1mm for cr R = 1.6/im), we find 
that they are good on the total space window (|z| up 
to 3000ct r ). As shown in Fig. 2(a), the extracted val- 
ues 7 fit (black dots) fairly agree with Eq. (2), except for 
low energy where the Born approximation breaks down. 
The values extracted in the same manner for each real- 
ization of the disordered potential are also shown (blue 
squares and red diamonds) . We find nonnegligible differ- 
ence between blue and red detunings (see Fig. 2), which 
can be ascribed to higher-order terms in the Born expan- 
sion [23]. Nevertheless, this difference is small and the 
strong increase of "/(E) appears for each realization in 
approximately the same region as predicted by Eq. (2). 
The parameters we used are relevant to current experi- 
ments as regards disorder [10, 12], observable space [10], 
and width of atom lasers [22]. It validates our proposal 
to probe the energy dependence of j(E). 

We now generalize the above results to higher dimen- 
sions (d > 1), for which the localization scenario is 
more involved. At intermediate distance (between the 
transport mean free path l B and the localization length 
L loc = I/7), where interference effects play a negligible 
role, the dynamics is diffusive with the diffusion constant 
D B (E) = (H/m)k E l B (E)/d [24]. For weak, isotropic dis- 
order [i.e. for c(k) = c(|«|)], one finds [25] 



m 2 Vlai 
(2Tr) d - 1 H 4 k 



3-d 



dSl d (l-cos0) c (2fc E CT R | sin(0/2)|) 

(6) 

with fid the hyperspherical angle in dimension d. On 
length scales larger than l B , interference effects induce 
AL, characterized by the Lyapunov exponent 7. The lat- 
ter can be calculated using the self-consistent theory of 
AL [26]. In 2D one finds, j(E) = l B 1 exp(-7rfc £ Z B /2). In 
3D, j(E) is the unique solution of [l — (7r/3)(fc E Z B ) 2 ] = 
7^ B x arctan(l/7Z B ), which exists only below the localiza- 
tion threshold (mobility edge), i.e. for k E l B < y3/V [25]- 
In both cases we can formally write j(E) — l^ 1 Fd(k E l B ) 
with Fa a decreasing function of k E l B , consistently with 
the scaling discussed in the introduction. It follows from 
Eq. (6) that, if c(/t) is as usual a constant or decreas- 
ing function of k, then l B (E) increases with E, and j(E) 
decreases when E increases. As for the ID case, this 
standard behavior can be changed by tailoring the disor- 
der correlations so that C(k) increases strongly enough 
in a certain window, and observed in the same way. 

In 2D, we propose to use an isotropic speckle potential 
created by a uniformly illuminated ring-shaped diffuser 
of inner radius r and outer radius R (see Inset of Fig. 3). 
For a thin enough ring (0.77R < r < R), we find that 
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Figure 3. (Color online) Lyapunov exponent versus energy 
in 2D speckle potentials created with a ring-shaped diffuscr 
of inner radius r and outer radius R (see Inset), for increasing 
values of r/R (from the upper line to the lower line on the 
left-hand-side of the figure) and \Vb\ = 0.25(fi 2 /mcr 2 .). 



j(E) is nonmonotonous with a marked local maximum, 
so that localization increases with energy in a given win- 
dow (see Fig. 3). For the parameters of Fig. 3, -f(E) 
peaks to about 5 x lO -4 ©-" 1 , with cr R = \ 1/2itR. For 
<7 a = 0.25/im, it corresponds to L loc ~ 500/im, which is 
within experimental reach [10]. Moreover, the width of 
the maximum is Afc ~ O.lu" 1 , which can be probed with 
the same atom laser as used in Fig. 2, the width of which 
is p„ — 0.0375/mt,^ 1 (note that a R is a factor 6.4 larger in 
the ID case above). 

In 3D, we consider an anisotropic scheme, easier to 
realize in practice. Then, Eq. (6) does not hold, but 
it can be anticipated that tailoring the correlations in 
anisotropic models of disorder can also invert the stan- 
dard behavior ofj(E). Consider the 3D speckle potential 
obtained by one Gaussian laser beam, of waist w (in both 
transverse directions) [same scheme as in Fig. 1(a)]. Due 
to the anisotropy of the disorder, the localization is de- 
scribed by an anisotropic Lyapunov tensor of eigenaxes 
x, y and z, and the Lyapunov exponent in all directions 
increases monotonically with E, as shown in Ref. [27]. 
If we now use two coherent parallel Gaussian beams of 
same waist w and separated by a distance A along Z 
(similarly as for the ID case), the interference between 
the two speckles create two bumps in C(k) at k ~ ±fcok z 
with ko = V2ttA/X 1 and k z the unit vector in the k z 
direction. These two bumps are expected to strongly en- 
hance the localization around the energy E cx T?k^j2m. 
Figure 4 shows the Lyapunov exponents found using the 
self-consistent theory of AL for anisotropic disorder [26- 
28]. As expected, the maximum is found for approxi- 
mately k E cx ko- For all the configurations of Fig. 4, 
7(2?) exhibits, below the right most mobility edge, a local 
maximum in each direction, hence realizing the desired 
effect. For the parameters of the right (purple) curve, 
j(E) vanishes (i.e. the localization length diverges) in a 
given energy window, thus opening a band of extended 
states inside the localized region, delimited by two new 
mobility edges. 

In summary, we have shown that the AL of noninter- 
acting quantum particles (matter waves) induced by a 
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Figure 4. (Color online) Lyapunov exponents along the 
transport eigenaxes (x: diamonds; y: dots; z: squares) and 
their geometric average (solid lines) versus energy in a 3D 
speckle potential created by two Gaussian beams (waist w 
and separation A) for various parameters (indicated in the 
figure). The arrows indicate the mobility edges. We recall 
that or = Xol/itw. 



correlated disorder in free space can increase with the 
particle energy in a given window. In contrast to other 
systems where this behavior is more common, e.g. elec- 
trons in crystal lattices or light waves in dielectric mate- 
rials, it is here purely due to appropriately tailored corre- 
lations of the disorder. We have proposed suitable meth- 
ods to tailor the correlations in optical disorder, which 
require moderate modifications of existing schemes in 
ID [10], 2D [12], and 3D [13, 14]. We have proposed a 
method to observe it in any dimension, which conversely 
differs from standard schemes used so far with ultracold 
atoms, and explicitly demonstrated its efficiency in the 
ID case. 

Let us finally discuss how the increase of j(E) with E 
can serve as a smoking-gun evidence of quantum local- 
ization of particles. For any experiment on localization, 
AL should be discriminated from other possible effects. 
For light waves for instance, it is necessary to distinguish 
it from absorption, which also produces exponential de- 
cay of the intensity. This can be done by analyzing the 
statistics of transmission [29, 30]. In contrast, ultracold 
atoms are not subjected to absorption, but they can be 
classically localized (trapped) in potential wells, below 
the localization threshold. Then, absence of diffusion 
and exponential decay of density profiles can hardly be 
viewed as indisputable proof of AL. For instance, classi- 
cal localization in some non-percolating media can lead 
to qualitatively similar effects, for instance in 2D speckle 
potentials [31]. For any model of disorder however, the 
classical localization length, defined as the average size 
of the classically- allowed patches [32], increases with the 
particle energy. Hence, the decrease of the quantum lo- 
calization length with the particle energy discussed in 
this work has no classical equivalent, and can be viewed 
as a smoking gun of quantum localization. This effect 
could be useful to demonstrate AL, in particular for 2D 
speckle potentials, which have a percolation threshold 
significantly higher than their 3D counterparts. ^From a 
practical point of view, it does not require accumulation 
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of many statistical data, in contrast to standard methods 
used for classical waves [29, 30]. 
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